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ABSTRACT 

The magnetic fields affect collapse of molecular cloud cores. Here, we consider 
a collapsing core with an axial magnetic field and investigate its effect on infall 
of matter and formation of accretion disk. For this purpose, the equations of 
motion of ions and neutral infalling particles are numerically solved to obtain 
the streamlines of trajectories. The results show that in non-steady state of 
ionization and ion-neutral coupling, which is not unexpected in the case of infall, 
the radius of accretion disk will be larger as a consequence of axial magnetic field. 

Subject headings: star: formation - ISM: clouds - ISM: evolution - magnetic 
fields 

1. Introduction 

A lot of observational information is now known about the structure of dense cores as 
progenitors of stars within the molecular clouds (e.g., Evans 2010). Theoretically, we expect 
that most of these cores will collapse into protostars, but the details of evolutions are less 
evident. Clearly, the evolution heavily depends upon the effects of local density, pressure, 
rotation, magnetic fields, the presence or absence of nearby stars and protostars, and upon 
the other physical phenomena. The most commonly used picture of low-mass star formation 
is inside-out collapse (Shu 1977) in which it begins at the center of a singular isothermal 
sphere and an expanding wave of infalling matter propagates outward with the speed of 
sound. Up to now, different models for core collapse and formation of protostars have been 
presented, in each the effects of some physical phenomena investigated (e.g., McKee and 
Ostriker 2007). 
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The velocity maps of molecular cloud cores, which are built by using the velocity mea- 
surements along the line of sight, show a small velocity gradient across each core. Measuring 
this global velocity gradient, with assumption that core has a uniform rotation and follow 
a rigid-body rotation law, can be used to deduce the core angular velocity (e.g., Goodman 
et al. 1993). The large size of initial collapsing core implies that even the modest initial 
rotational velocities will cause the infalling matters to land first on a rotationally supported 
disk rather than a pressure-supported proto-star (e.g., Hartmann 2009, Nejad-Asghar 2011). 
In the simplest analysis of rotating collapse, we assume that pressure forces are negligible 
and so the problem can be analyzed by using ballistic trajectories. The results for collapse of 
a spherically symmetric cloud in uniform (solid-body) rotation were initially worked out by 
Ulrich (1976), with subsequent extension to disk formation by Cassen and Moosman (1981), 
and to collapsing singular isothermal sphere by Terebey, Shu and Cassen (1984). Mendoza, 
Tejeda and Nagel (2009, hereafter MTN) have recently generalized this idea by construction 
a steady analytic accretion flow model for a finite rotating gas cloud. They show that the 
streamlines and density profiles deviate considerably from ones calculated by Ulrich, and for 
all relevant astrophysical situations, the assumption of a finite cloud radius needs to be used. 

The observations indicate that magnetic energy in the molecular clouds is comparable 
to the gravitational energy (e.g., Crutcher 1999). Moreover, the magnetic fields are also 
theoretically believed to play an important role in gravitational collapse of the molecular 
cloud cores. They provide pressure support against the gravity and carries away angular 
momentum prior to and during the collapse of cores to form accretion disks, jets and pro- 
tostars (e.g., Machida 2010). Although, in a molecular cloud core, the spatial configuration 
of magnetic field lines is not simple (e.g., Whittet 2005), but polarimetry observations of 
young stellar objects suggest that circumstellar thin disks around the young stars have ap- 
proximately aligned perpendicular to the magnetic fields (e.g., Pereyra et al. 2009). Here, 
we consider a simple initial configuration in which the magnetic fields are assumed to be 
parallel with rotational axis of core. Since the molecular cloud cores are lightly- ionized, the 
ambipolar diffusion in which the magnetic field is frozen into the charged species and drifts 
along with them through neutrals, is an important mechanism (e.g., Adams 2009). Thus, we 
consider the effect of magnetic fields directly on charged particles, while the neutral species 
feel them indirectly via the collisions with ions. 

In this way, as an extension to the work of MTN, we investigate the effect of axial 
magnetic field on streamlines of infalling ions and neutral particles. For this purpose, formu- 
lation of problem and equations of motion are given in section 2. In section 3, trajectories of 
infalling ions and neutral particles are investigated and effect of the magnetic fields on the 
accretion disk radius is studied. Finally, section 4 is allocated to summary and conclusions. 
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2. Formulation of problem 



The molecular cloud cores are lightly-ionized gases with small fraction of ions with 
density rii and electrons with density n e ~ rij. If the volumetric rate of electron- ions re- 
combination took place in gas phase, it would be proportional to n e rii oc nj. In a simple 
steady state approximation, the volumetric rate of electron-ions recombination oc nf would 
be equal to the volumetric rate of neutrals ionization via cosmic rays oc n n . Although, the 
actual ionization balance in the molecular clouds is more complicated, Elmegreen (1979) and 

1 /9 

Umebayashi and Nakano (1980) showed that the relation n-i oc n n is a good approximation. 
Here, we suppose for fiducial purpose that pi = eepl/ 2 , where e = 9.5 x 10 _15 m _3 / 2 kg 1 / 2 and 
e is a dimensionless free parameter which represents the deviations from calculations in the 
steady state approximation. 

In the magnetized cloud cores, the mean velocity v n of neutrals will not generally equal 
to the mean velocity of ions and v e of electrons. Although, the positive ions and electrons 
feel the magnetic forces in opposite direction, but the induced electric field will generally 
cause they move in ambipolar motion, so that in the time-scales that are considered here, 
their mean velocities are approximately the same (v e ~ Vj). Since electrons carry much less 
momentum than ions, we neglect the dynamic of electrons. In this way, resistance of rela- 
tive drift between ions and neutrals will be a drag acceleration arises from mutual collisions 
between them. The drag acceleration exerted on neutrals by ions is 77Pi(v; — v„); the drag 
acceleration exerted on ions by neutrals is — 77p„(vj — v n ), where 7 = 3.5 x 10 10 m 3 kg _1 s _1 
(Draine, Roberge and Dalgarno 1983) is the drag coefficient and 7 is a dimensionless free 
parameter which represents the deviations from calculations in the steady state approxima- 
tion. 

The equations of motion for ions and neutral particles, in gravitational field of a central 
mass, are respectively as follows 



where M is the proto-stellar mass at the origin, and = N\e and rrii = N 2 m p are ion charge 
and mass, respectively, as an integer multiple of the electron charge e and the proton mass 
m p . Here, we assume that N 2 ~ 2N X . 

We measure the length and time in units of O.Olpc and 10 5 yr, respectively, so that 
velocity unit is 0.1 km s" 1 . Mass unit is appointed to 1M©, thus the gravitational constant 
is G = 45. In this manner, we have 7 = 7.8 x 10 9 and e = 3.5 x 10 -8 . Finally, the unit 




GM 





n 
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of magnetic field is chosen equal to 10 nT = 100 [iG. Using these units, the equations of 
motion ([T|) and ([2]) can be rewritten in a dimensionless form as follows 

M 

r- = -45— n + 0.5Sr 4 x B - 3.9 x lO^n^ - r„), (3) 

' i 

i n = -45^r n + 1.9e7ny 2 (r 4 - r n ), (4) 

where B is a unit vector in the direction of magnetic field, and n n = n ra /(10 10 m -3 ) and B 
are two dimensionless free parameters. 

Choosing the axial magnetic field B = k, the equations of motion in Cartesian coordinate 

are 

oti = —45M sin 8i cos (pi— + 0.5Byi — 3.9 x 10 5; yn n (xi — x n ), (5) 

^ i 

y\ = -45M sin ^ cos 0— ~~ 0.5-Bx, - 3.9 x 10 s jn n (yi - y n ), (6) 

i 

Zi = — 45Msin6 l j cosfa-^ ~~ 3-9 x 10 5 7n n (ij — z n ), (7) 



x„ = -45Msin6 l n cos0 n — + OeynJ/ 2 ^ - x n ), (8) 



n 



jj n = -45M sin # n cos + lMju]/ 2 ^ - y n ), (9) 



£ n = -45Msin0 n cos0 n -r + 1.9e7n^ /2 (i, - i„), (10) 



where 9i, 9 n , fa, and n are the polar and azimuthal angles of ions and neutral particles, 
respectively, and r» = \/ xf + yf + z 2 and r n = -\/:r 2 + + z\ are their radial distances. 
There are four free parameters B, e, 7, and n n which represent the strength of magnetic 
field, ionization fraction, drag coefficient, and neutral density, respectively. 



3. Trajectories of infalling matters 

As mentioned before, there are many observations that show velocity gradients in the 
maps of molecular cloud cores. These velocity gradients can be used to deduce the rotation 
of cores, but evaluation of their angular velocity Qo is considerably complicated because 
the velocity fields usually exhibit complex supersonic motions. Here, we choose the fiducial 
value of Qq = 3 x 10 _13 rad s _1 , which in the time unit 10 5 yr is f2 ~ 1. The effects of initial 
angular momentum in core collapse and formation of accretion disks are now approximately 
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well known (e.g., Hartmann 2009). Here, we turn our attention to the effects of magnetic 
fields and drag force on infalling matters. If we assume that 7 = and B = 0, equations 
©-(El) lead to the trajectories which are contained on a plane and given by a conic section 
(i.e., the results obtained by MTN). 

Choosing the initial conditions 
r n = n = r r 

= r sin 6 Q cos0o* + r o sin 9 sin o j + r o cos 9 k, (11) 

v„ = Vj = -v 0r f + v 0<t) (f) = -(vor cos O + r sin O ) sin 6 i 

— (vor sin 0o — r o cos <f>o) sin 9 j — v 0r cos 9 k, (12) 

for position and velocity of infalling particles, the equations of motion (|5|)-( IT0|) can be 
solved by numerical methods such as Runge-Kutta. Approximately, all cores in the maps 
of molecular clouds seem apparently to be elongated rather than spherical (e.g., Curtis and 
Richer 2010). Since determining the exact three-dimensional shape of a core from apparent 
observations of the-plane-of-sky is impossible, statistical techniques have to be applied; as 
the results, some works indicate a preference for prolate cores (e.g., Gammie et al. 2003, Li 
et al. 2004), and some others favour oblate shapes (e.g., Curry 2002, Jones and Basu 2002). 
Theoretically, the prolate elongation of the cores may be inferred as a remnant of their origin 
in filaments and turbulent flows (e.g., Hartmann 2002), while the oblate shape of collapsing 
cores may be justified from the effects of magnetic fields and rotational motion (e.g., Nejad- 
Asghar 2010). The streamlines of infalling particles, in absence of magnetic field (B = 0, 
7 = 0) and initial radial velocity (t>o = 0), for initial spherical core (r = 1), oblate core 
(r = 1 / a/ sin 2 #0 + 4 cos 2 6 ), and prolate core (r = 1/ a/4 sin 2 8 + cos 2 6 ) are shown in 
Figffl 

The numerical plots of accretion disk radius for neutral particles by substitution 6q — > 
7r/2, for a wide variety of initial parameters, are shown in Fig. [2J The numerical results 
show that the accretion disk radius is approximately proportional to the inverse of initial 
radial velocity vo r as was given by the analytical ballistic models of MTN. Since we want to 
investigate the effect of magnetic fields on trajectories of infalling matters, without losing 
the generality of problem, we consider an initial spherical core with assumption vq t = 0. 
Here we turn our attention to only one streamline with r = 1 and 6 = 7r/4. The effects of 
axial magnetic field on this streamline are shown in FigfS] for three different values of free 
parameter B with a fixed value of three parameters 7, e and n n equal to 10 -5 , 10 5 , and 
1, respectively. As a result, we see that the axial magnetic field causes the infalling ions 
land at a larger equatorial radius, and collisions of ions with neutral particles lead to the 
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same behavior for infalling neutral matters. In FigJU the equatorial radii of ion and neutral 
infalling particles (from 9$ = 7r/4) are depicted versus 7 and e at a fixed value of magnetic 
field. We see that increasing of drag coefficient 7 leads to coupling of ions and neutrals so 
that their equatorial radius reach to each other asymptotically. The same behavior can be 
seen for the effect of ionization degree e so that increasing of this parameter asymptotically 
causes to approximately the same equatorial radii of ions and neutral infalling particles. The 
third free parameter is neutral density h n which its effect is the same as 7 but its effect on 
neutral particles is less. 

The density of neutral particles can be evaluated by assuming that the mass infall rate, 
M, is steady. Assuming that the cloud at r is nearly spherical, the mass flow in a flow tube 
spanned by A9 is 

. ■ 2ixr 2 sin 9 MA9 1,. „ A „ 

AM = = -Msin# A# . 13 

Anr 2 2 

The density at (r, 9) can be found by following the streamlines corresponding to 9q and 
9 + A9 , 

p(r , 9) = ^ = _JL_^^±. (14 ) 

27rr 2 sin6A6 l |tv| 47rr 2 |tv| sm.9 A9 

In order to obtain the iso-density contours, we divide < #0 < it /2 to N equal segments 
so that A9 = tt/2N. For two streamlines in the borders of each segment, values of sin#o, 
sin# and v r are evaluated in the mid points. In this way, the density at different radii 
and polar angles can be found from ([Hj), and the iso-density curves can be depicted as 
typically are shown in Fig. In this figure, the mass accretion rate is assumed to be 
M = 4tt x 10~ 7 Mq yr _1 . These iso-density contours show that the axial magnetic fields 
cause the infalling matters land at a larger equatorial radius. 



4. Summary and conclusions 

Nowadays, we have accepted that the molecular cloud cores can be gravitationally un- 
stable under certain conditions, and they can collapse. Rotational motion of the cores leads 
to land the infalling matters on an accretion disk around the proto-star. The streamlines of 
infalling matters for initial spherical, oblate and prolate cores are depicted in FigJTJ Obvi- 
ously, the centrifugal force causes to increase the initial dimension of accretion disk in oblate 
cores, while it decreases the disk radius in prolate ones. 

Next, we considered an axil magnetic field in the collapsing core. For investigation the 
effect of this magnetic field, we wrote separately the equations of motion of neutral and ion 
infalling particles. In general form, these equations are presented with four free parameters 
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B, e, 7, and h n which represent the strength of magnetic field, ionization fraction, drag 
coefficient, and neutral density, respectively. The accretion disk radius for different values 
of the relevant input parameters are depicted in FigfSJ Since our goal is to study the effect 
of magnetic field and drag coefficient on infalling matters and disk formation, we turned our 
attention to only one streamline at 9q = 7r/4; the results are shown in Figj3j 

In a steady state of molecular gas, the ions are well-coupled to neutral particles (i.e., 
e = l,7 = l) so that their dynamics are approximately the same. Fig JH shows that increasing 
of drag coefficient and ionization degree lead to the same behavior for ions and neutral 
infalling particles. In infalling case, the frequency of collisions between ions and neutral 
particles eventually becomes so low that the balance cannot be maintained, thus, the values 
of free parameters e and 7 deviate from the steady sate. In this non-equilibrium case, ions 
are directly affected by magnetic field so that they land in a larger radii at equatorial plane. 
In the low-coupling (7 << 1) and high-ionization degree (e >> 1), the dynamics of ions 
can affect the motions of neutral particles so that they also land in a larger radii to form 
the larger accretion disk as shown in Figj3j The same behavior for increasing of accretion 
disk radius via the effect of magnetic field can be deduced from iso-density contours that are 
depicted in Fig. [5J 
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Fig. 1. — The streamlines of ions and neutral particles in absence of magnetic field and initial 
radial velocity (B — 0, 7 = 0, v 0r = 0), for (a) initial spherical core (r = 1), (b) initial oblate 



core (r = 1/ysin 6 + 4 cos 2 6*o), and (c) initial prolate core (r = l/v4sin 6*o + cos 2 ^o)- 
In each figure, the curves from top to bottom have the initial azimuthal angle equal to 
9 = 15, 30, 45, 60, and 80°. 
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Fig. 2. — Accretion disk radius by the substitution 6 — > it /2 for various initial free param- 
eters with h = 1. 
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(b) 

Fig. 3. — The streamlines of (a) neutral and (b) ion flows for different values of free parameter 
B = 1 (solid), B = 5 (dash), and B = 10 (dot), at 6 = vr/4 with 7 = 10" 5 , e = 10 5 and 
n n = 1. 
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Fig. 4. — The equatorial radius of infalling ions and neutral particles which start at 6q = 7r/4, 
versus (a) logarithm of free parameter 7 with e = 10 5 , and (b) logarithm of free parameter 
e with 7 = 10~ 5 . In these figures the free parameters B and n n are chosen equal to 5 and 1, 
respectively. 
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Fig. 5. — The typical iso-density contours for p « 0.01 (top panel), and p ~ 0.1 (bottom 
panel), for different values of free parameter B = 1 (solid), B = 5 (dash), and B = 10 (dot), 
with 7 = 10" 5 , e = 10 5 and h n = 1. 



